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Abstract 

We obtain a complete set of one-loop RGE's for a set of combinations of neutrino parameters 
for the case of two-fold degenerate hierarchical three-neutrino models. The requirement of consis- 



tency of exact solutions to these RGE's with the two-fold degeneracy yields conditions which have 
previously been obtained perturbatively /numerically. These conditions, in the limit \Ueuz \ = 0, are 
■ shown to lead to a strong cancellation in the matrix element of neutrinoless double beta decay. 
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It has been a puzzle that the mixing in the leptonic sector is so large while the mixing 
in the quark sector is so small. Many attemps have been made to explain this fact. One 
possible scenario is to utilize the flavour symmetry combined with GUT symmetry at some 
high energy scale A [|I|. Viable symmetries are those giving rise to large mixing in the 
lepton mass matrices. Most models of this kind suffer from fine-tunning and the difficulty 
of constructing a viable superpotential in the flavour symmetry sector that gives rise to the 
required vacua. An alternative to this scenario is the idea of infrared fixed point (IRFP) P|. 
Contrary to the idea of flavour combined with GUT symmetry, in the IRFP scenario, the low 
energy physics is governed by the low energy dynamics, namely, the renormalization group 
equations below the scale A. Physics above the scale A plays no role in the predictions at 
low energies. Therefore, if there exists any IRFP which leads to viable phenomenology, one 
does not have to deal with the fine-tunning problem and the difficulty of finding the correct 
vacua. The focus of our attention in this note is the hierarchical three neutrino models with 
two-fold denegeracy, and the implications of the exact solutions to one-loop renormalization 
equations (RGB's) to these models. 

In the flavour basis where the charged leptons are diagonal, the neutrino flavour eigen- 
states and mass eigenstates are related by >= Uail^i >, where a and i are the flavour 
and mass eigenstate indices respectively. The mass matrix can be diagonalized as follows 

U'^niyU = diag{mi, 1712, rris) (1) 

We adopt the usual parametrization for the leptonic mixing matrix U 

U = diag{e'^%e'^^,e'^^) ■ V ■ diag{e-''^''^ , e'^^'/^, 1) (2) 



V 



^ C2C3 C2S3 S2e-'^^ 

-C1S3 - SiS2C2,e^^ C1C3 - SiS2C3e*'^ S1C2 

y S1S3 - CiS2C3e*'^ -S1C3 - CiS2S3e*'' C1C2 y 



(3) 



where Cj = cos6'j and Sj = sin^j and < < 7r/2. The 5e.^i,r are three unphysical phases 
which can be absorbed by phase redefinition of the neutrino flavour eigenstates. There are 
three physical phases: 5 is the universal phase (analog of the phase in the CKM matrix), 
and and (j) are the Majorana phases. By properly choosing the phases and (j) all 
three mass eigenvalues rrii can be made positive. We therefore assume, without loss of 



generality, that (mi,m2,m3) are positive. If any of these three phases is not zero or not vr, 
CP violation in the lepton sector is implied. Note that in the limit 62 = 0, 9i is identified 
as the atmospheric mixing angle, Oatm-, and 9^ is identified as the solar mixing angle, ^0. In 
general, the mixing matrix elements are related to the physical observables, the atmospheric 
and solar mixing angles, hy svo? 26atm = 4|f/^3p(l — |f/^3p), and sin^ 26*0 = 4|[/e2p(l — |t/e2p)- 
Recent results indicate that for atmospheric neutrino oscillations, Am^^^ = 3.1x 10~^ey^, 
sin^ 29atm = 0.972 [^]; for solar neutrino anomaly problem, there exists four solutions: (i)VO: 
Am| = 8.0 X 10-"e1/^ sin^ 2^0 = 0.75, (ii)LOW: Am| = 7.9 x lO^^eV^, sin^ 2^0 = 0.96, 
(iii)LAMSW: Am| = 1.8 x IQ-^eV^, sin^ 2^0 = 0.76, (iv)SAMSW: Am| = 5.4 x lO-^eV"^ 
sin^ 26'0 = 6.0 x 10^'^ |Q ; and the matrix element \Ue3\ = sin 62 is constrained by the CHOOZ 
experiment to be jt/esl < 0.16 0. 

The observed relation Am^^^ = \m^ — mW ^ Amg = \m\ — m\\ in the two- 
fold degenerate, hierarchical model implies ^ m2 — rrii, and this in turn implies 
V21 > V32 ^ V31 ~ 1 with Vjj = {rrii + mj){mi - mj)^^. 

We assume that the neutrino masses are generated by a dimension-5 effiective Majorana 
mass operator in the MSSM 



The neutrino mass matrix {'m,y)ij is related to kij by m,y = kijv'^/2, where v"^ = v"^ + 
v1 = (246ey)^ is the squared vacuum expectation value of the SM Higgs. The effective 
dimension-5 operator is generated by some mechanism at the high energy scale A. The 
seesaw mechanism is the most common way to generate this operator. Since we are only 
interested in physics below the scale A, we will start with the effective Lagrangian Eq. @ 
without specifing the origin of this effective operator. 

The general one-loop RGE of the effective left-handed Majorana neutrino mass operator 
is given by 



£ D -kij{HuLi){HuLj) + h.c. 



(4) 



dm, 



(5) 



dt 
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where t = In /x. In the MSSM, P and /t„ are given by, 



1 y^y 1 

SzTT^ cos^ p 327r^ cos"' p 
- -[7^1+6^2-6- 



- Tfr^[|^i+ 6^^2-6^] (7) 
loTT^ 5 sm p 

where (/f = l^fy is the U{1) gauge couphng constant, F„ and Yg are the 3x3 Yukawa couphng 
matrices for the up-quarks and charged leptons respectively, and ht and h^- are the SM t- 
and r- Yukawa couphngs. Since gives rise to an overall rescaling of the mass matrix, it 
has no effects on the running of the mixing matrix U. Eq. can be solved analytically by 
integrating out its right-hand side [§. Note that at one-loop level, since the evolutions of 
the gauge coupling constants gi,2if) and of the diagonal Yukawa couphngs ht^rif) are known, 
it is indeed possible to carry out the integrations on the right-hand side without making 
any further assumptions. However, the diagonalization procedure of the resulting 3x3 
complex symmetric matrix, m^{t), is very complicated. It is thus hard to infer analytically 
the behaviours of the physical observables, the mixing angles and phases. An alternative to 
this "run-and-diagonalize" procedure is the "diagonalize-and-run" procedure. It is convenient 
to work with the RGE's of mass eigenvalues and the diagonalization matrix, given by ^ 

= -2miPii - miRe{Ku} (8) 

(9) 

where 

Tii = iQii (10) 



= {— ^){{ml + m])Pij+2mimjP*^} + iQij 



— TTij 



Vi,Re{Pij} + tVr^'lmiPij} + zQ^j (11) 



Here P and Q are defined as 



P = ^U\P + P^)U, Q = ^U\P-P^)U (12) 
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Eq- ®, (I), and (H) have been studied before [|, |T^, 0, [T^, |T3|, [T^, but the analyses 



have been done either numerically or perturbatively. (Exact solutions to the RGE's in 
a two flavour case have been investigated recently in [Q). Due to the large interfamily 
hierarchy in the charged lepton sector, we keep only the r- Yukawa coupling. We will further 
assume that hr does not evolve throughout the entire range of the RG running. This is a 
valid assumption for the hierarchical case with two-fold degeneracy, as V21 is very large. 
Under these assumptions, the above quantities are given in terms of the masses and the 
diagonalization matrix elements as 



4 = -^f^3.f^3„ Q = (13) 



The evolutions of Vj^, {Pa — Pjj) and Re{Pij) can be derived from Eq. (H) and (^. For 
(i,j) = (2,1), with V21 ^ V31 — V32 — 1, the RGE's for these three functions form a 
complete set of coupled differential equations as follows: 

^^^^ =V^,(P22-Pn) (14a) 



dt 

^ -4V2i[i?e(P2i)]' (14b) 
V2i(P22-Pii)i?e(P2i) (14c) 



diP22 - Pu) _ ,,2 

dt 

dRe{P2i) 



dt 

The exact solutions to these coupled differential equations are given by 

V2i(t) = a^Z{t)~^'^ (15a) 
(P22(t) - Pii(t)) = {hi + 4cg(l - Z{t)-^) f'' (15b) 
Re{P2i{t)) = coZ{t)-^'^ (15c) 

where 

Z{t) = 1 - 2aobot + alibi + 4cl)t^ (16) 
and ao, &o and cq are the initial values at the high energy scale A: 

ao = V2i(0); 6o = (H2(0)-Pii(0)) 

Co = Pe(P2i(0)). (17) 

The behaviours of these three functions are shown in Fig. (|l])-(^. Note that V2i(t) and 
Re{P2i{t)) flow to zero, while (P22(i) — -Pii(i)) flows to a constant value of (69 + 4co)^/^ in 



the infrared. This set of parameters, (V2i(t*), i?e(P2i(t*)), (H2(t*) - Pnit*))) = (0, 0, {b^ + 
4co)^^^) is an infrared stable fixed point; however, it is unreahstic. The function V2i(t) 
decreases to 0(1) very fast as the energy scale goes down, for any non-vanishing Bq and Cq, 
however small they are. This is phenomenologically unacceptable. In addition, it contradicts 
with the assumption V21 ^ V32,V3i we made in order to arrive at Eq. ([l^)- For the 
consistency of the calculations, we thus require the following two conditions at the initial 
high energy scale A: 

bo = (P22(0)-Pii(0)) = (18) 
Co = /2e(P2i(0)) = (19) 

We emphasize that these conditions have been obtained by demanding that the exact so- 
lutions to the above RGB's Eq. (|1^) be consistent with V21 ^1. It is to be noted that 
these conditions have been obtained before numerically and perturbatively ||^, |rT|. When 
these conditions are satisfied, all three equations in Eq. (|1^ do not evolve. The first relation, 
Eq. (p!8|), gives rise to IV32P = IV31P which translates into 

22 2 

''^''^-^i -tan 2^3 -cos 5 (20) 



sin 29i ■ S2 

The second relation, Eq. (|19D, gives rise to 



sin 26*3 {cfsl — sf) cos 26^3 
sin 6 ■ sin 26i ■ S2 tan 6 

Combining these two relations, we obtain a very simple relation among 6^ and three CP 
violating phases S, 0, 0': 

cos 2^3 = - l—r- (22) 

tan 5 tan(^) 

We have studied the RGE's involving various functions, Vjj, {Pa — Pjj) and Re{Pij), for the 
case = (3, 1) and (3,2). Upon imposing the above consistency conditions Eq. ([18| ) and 
(pj|), we deduce that the functions Pn, P22, P33, Re{P3i) and Re{P32) do not run. These 
results cannot be tested experimentally at present. 

Now we discuss the implications of Eq. (18) and (|19D in the limit 62 = (Recently, it 



has been pointed out that this could be a consequence of the so-called 2-3 symmetry 
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FIG. 1: The function V2i(i)- Initial values at A are (ao, bQ, cq) = (1000, 1, 1) 



They imply 



cos 



)SiC3S3 = 0, 



2/ 2 







(23) 



Since the atmospheric angle 6'i = 7r/4 is non- vanishing, these two relations can be satisfied 
simultaneously only if (i) the solar mixing angle is maximal, i.e. 6^ = 7r/4, and (ii) the 
Majorana phase difference (0 — (p') = n. The phases (p and (p' occur in the matrix element 
(Mee) for the neutrinoless double beta decay: 



j=l,2,3 



(24) 



where index i denotes the mass eigenstates. Currently, the most stringent bound is given 



by i? = 0.2eV [T^. In the limit 62 = with nearly degenerate mi ~ m2, M^e becomes 



Mee ^ \mi4{e-"^ci + e-''^ si) 



(25) 



It is obvious that when (0 — (f)') = tc and ^3 = tc/A the r.h.s. of Eq. ( p5D is exactly zero. 
Thus we conclude that neutrinoless double beta decay is very highly suppressed. 

It is interesting to speculate the reason(s) for consistency conditions of Eq. ([T8| ) and (|19|). 
It could be due to the existence of a symmetry at a high energy scale A. The other possibility 
is that these two relations are the fixed point relations of the RGE's for new physics above 
the scale A. 
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FIG. 2: The function (P22(*) - Puit))- Initial values at A are (ao, 6o, cq) = (1000, 1, 1) 
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FIG. 3: The function Re{P2i{t)). Initial values at A are (ao,6o,co) = (1000, 1, 1) 
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